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ABSTRACT 


The  effects  on  the  pressure  and  heat-transfer  distribu¬ 
tions  of  injecting  helium  and  argon  into  an  air  laminar 
boundary  layer  over  a  spherically  blunted  7.5-deg  half-angle 
cone  at  Mqq  =  13.2  and  Reoo/in.  =  1520  were  considered.  Calcu¬ 
lations  are  presented  for  a  binary  perfect  gas  boundary  layer 
including  transverse  curvature  and  displacement  effects,  an 
ionizing  air  nonequilibrium  boundary  layer  without  higher- 
order  effects,  and  a  nonequilibrium  stagnation  point  thin 
viscous  shock  layer  model.  The  effects  of  inviscid  chemical 
equilibrium  and  nonequilibrium  external  conditions  and  wall 
catalyticity  were  investigated.  Neglecting  all  higher-order 
effects,  comparison  of  binary  gas  and  ionizing  air  boundary- 
layer  results  showed  reasonable  agreement.  Strong  effects  of 
inviscid  chemical  nonequilibrium  on  stagnation  point  heat 
transfer  were  found.  Arguments  are  given  for  using  inviscid 
nonequilibrium  stagnation  conditions  rather  than  the  usual 
chemical  equilibrium  conditions  for  nonequilibrium  boundary- 
layer  calculations. 
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NOMENCLATURE 

Unless  otherwise  noted,  all  lengths  are  nondimensionalized 
with  respect  to  the  nose  radius  (see  Fig.  1,  Appendix  I,  for 
cone  geometry) 

Cj  Skin-friction  coefficient 

00 

Specie  mass  fraction 
Cp  Constant  pressure  specific  heat 

Binary  diffusion  coefficient 
ECW  Equilibrium  catalytic  wall 

FCW  Fully  catalytic  wall 

f'  Velocity  ratio,  u/ue 

f"  Velocity  gradient,  d(u/ue)/dn 

g  Stagnation  enthalpy  ratio 

H  Stagnation  enthalpy 

h  Static  enthalpy 

j  Diffusive  mass  flux 

k  Thermal  conductivity;  reaction  rate 

Le  Lewis  number 

M  Mach  number 

m  Mass  flow  rate 

NCW  Noncatalytic  wall 

NEQ  Nonequilibrium 

Pr  Prandtl  number 

p  Static  pressure 

p '  Normal  shock  stagnation  pressure 
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q  Heat  flux 

R  Gas  constant 

Re  Reynolds  number 

Rr  Nose  radius,  1  in. 

r  Radius 

• 

5  Net  rate  of  mass  production  due  to  surface  reactions 

s  Surface  distance  from  the  forward  stagnation  point 

T  Temperature 

TVC  Transverse  curvature 

U  Free-stream  velocity 

00 

u  Tangential  velocity  component 

v  Normal  velocity  component 

x  Surface  distance 

y  Distance  normal  to  the  surface 

Z  Air  specie  mass  fraction 

a  Local  body  angle 

7,7  Ratio  of  specific  heats 

A  Shock  layer  thickness 

6  Boundary-laye^  thickness 

6*  Boundary- layer  displacement  thickness 

6*^  Boundary-layer  displacement  thickness  with  mass  transfer 

e  Density  ratio,  p^/Pg 

•q  Transformed  distance  normal  to  the  surface 

Q  Temperature  ratio,  T/Tg 

m  Dynamic  viscosity 
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p  Density 

t  Shear  stress 

SUBSCRIPTS 

e  Edge  of  boundary  layer 

eff  At  the  effective  body  surface 

f  Outer  (air)  species 

i  Injected  species 

o  Stagnation 

s  Normal  shock  conditions 

TVSL  Thin  viscous  shock  layer 
w  Wall 

oo  Free-stream  conditions 
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SECTION  I 
INTRODUCTION 


The  study  at  the  Arnold  Engineering  Development  Center 
(AEDC) ,  von  Karman  Gas  Dynamics  Facility  (VKF)  of  effects  of 
mass  transfer  into  laminar  boundary  layers  at  hypersonic  wind 
tunnel  conditions  began  several  years  ago  with  the  work  of 
Whitfield  and  Griffith  (Ref.  1).  They  investigated  the  ef¬ 
fects  of  tangential  injection  of  helium  and  nitrogen  on  the 
zero-lift  drag  of  a  9-deg  half-angle  cone  at  =  10  and  20. 
The  effects  of  gas  injection  normal  to  the  surface  were  in¬ 
vestigated  by  Lewis,  Marchand,  and  Little  (Ref.  2)  for  sharp 
9-deg  half-angle  cones  at  Mw  -  9  and  a  5-deg  half-angle  cone 
at  =  3.93  and  5.64.  In  this  latter  work,  the  analysis  of 
mass  transfer  effects  on  zero-lift  drag  was  restricted  to 
normal  injection  of  air  or  nitrogen  into  an  air  or  nitrogen 
boundary  layer  based  on  the  method  of  Clutter  and  Smith  (Ref . 
3).  The  effects  of  injecting  helium,  argon,  and  carbon 
dioxide  into  air  or  nitrogen  sharp  cone  laminar  boundary 
layers  were  investigated  by  Mayne,  Gilley,  and  Lewis  (Ref.  4) 
using  the  boundary-layer  theory  of  Jaffe,  Lind,  and  Smith 
(Ref.  5).  In  Ref.  5,  a  binary  mixture  of  perfect  (nonreact¬ 
ing)  gases  was  considered,  and  the  effects  of  transverse 
curvature  (TVC)  were  correctly  included  whereas  the  treat¬ 
ment  in  Ref.  3  was  only  approximate.  (See  Ref.  6  for  a  dis¬ 
cussion  of  the  effects  of  the  TVC  approximation.) 

The  effects  of  mass  transfer  into  a  nonreacting  air 
boundary  layer  were  considered  by  Lewis,  et  al .  (Ref.  7) 
based  on  the  methods  of  Jaffe,  Lind,  and  Smith.  The  effects 
of  nonequilibrium  chemistry  on  the  inviscid  flow  field,  and 
the  no-injection  ionizing  air  laminar  boundary  layer  were 
also  considered  in  Ref.  7.  The  nonequilibrium  laminar 
boundary-layer  theory  used  in  this  most  recent  study  was  an 
extension  of  the  method  developed  and  described  by  Blottner 
(Ref.  8). 

The  theoretical  models  and  numerical  techniques  used  in 
the  previous  studies  have  been  extended  and  modified  for  use 
in  calculations  to  be  discussed  in  this  report.  The  theoret¬ 
ical  models  will  be  described  in  the  following  section,  and 
then  the  conditions  used  in  the  calculations  will  be  given. 

A  presentation  and  discussion  of  the  results  of  the  calcula¬ 
tions  will  be  given  followed  by  the  concluding  remarks. 
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SECTION  II 

THEORETICAL  MODELS 


The  inviscid  gas  models  used  to  describe  the  outer  flow 
field  are  presented  and  discussed  below.  Next,  the  binary 
perfect  gas  and  multicomponent  ionizing  air  boundary- layer 
models  are  described.  Finally,  a  thin  shock  layer  model  is 
described  for  the  ionizing  air  stagnation  region  with  injec¬ 
tion  of  either  helium  or  argon. 


2.1  INVISCID  FLOW  FIELD  MODELS 

The  boundary-layer  outer  edge  properties  were  obtained 
from  either  ideal  gas  (7  =  1.4)  blunt  body  and  characteristics 
solutions  or  from  a  nonequilibrium  blunt  body  solution  and  a 
streamtube  expansion  using  an  ideal  gas  (7  =  1.4)  blunt  body 
and  characteristics  solution  surface  pressure  distribution. 

The  binary  perfect  gas  laminar  boundary- layer  calculations 
used  the  ideal  gas  blunt  body  and  characteristics  solutions 
and  the  ionizing  air  boundary- layer  calculations  used  the 
nonequilibrium  blunt  body  and  streamtube  expansion  data. 


2.1.1  Ideal  Gas  Model 

The  ideal  gas  (7  =  1.4)  blunt  body  and  characteristics 
solutions  for  the  geometric,  spherically  blunted  cone  were 
obtained  by  the  method  developed  by  Inouye,  Rakich,  and 
Lomax  (Ref .  9) .  The  blunt  body  region  was  solved  by  the  in¬ 
verse  method  described  by  Lomax  and  Inouye  (Ref .  10)  in  which 
the  bow  shock  wave  location  and  shape  were  iterated  until  the 
zero  streamline  satisfactorily  matched  the  blunt  body  (sphere) 
geometry.  The  numerical  scheme  is  very  stable,  and  blunt 
body  solutions  were  obtained  in  about  three  minutes  and  the 
characteristics  solutions  required  about  15  minutes  on  a 
CDC  1604B  digital  computer  with  an  on-line  printer. 

The  blunt  body  and  characteristics  solutions  over  the 
geometric  body  perturbed  by  the  displacement  thickness  (in¬ 
cluding  mass  transfer  effects)  were  obtained  by  the  method 
described  by  Marchand,  Lewis,  and  Davis  (Ref.  11).  The 
method,  an  extension  of  the  method  of  Inouye,  Rakich,  and 
Lomax  (Ref.  9),  provides  a  blunt  body  solution  over  the 
sphere  perturbed  by  the  displacement  thickness  (an  ellipsoid) 
and  a  least-squares  curve  fit  of  the  downstream  "effective 
body"  defined  by  reff  =  rw  +  6*  cos  a.  The  method  has  been 
found  to  yield  accurate  fits  of  the  effective  body  geometry 
and  smooth  numerical  results  over  the  entire  perturbed  body. 
Some  details  regarding  the  procedure  used  and  some  results 
from  no- injection  calculations  are  given  in  Ref.  11. 
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2.1.2  Inviscid  Nonequi librium  Model 


Only  the  nonperturbed  inviscid  nonequilibrium  flow  field 
solution  was  used  in  this  study.  The  theoretical  model  used 
was  the  inverse  blunt  body  method  of  Curtis  and  Strom  (Ref. 

12)  to  obtain  the  nonequilibrium  inviscid  stagnation  point 
boundary- layer  edge  conditions,  and  the  nonequilibrium  stream- 
tube  program  of  Lordi,  Mates,  and  Moselle  (Ref.  13)  was  used 
to  obtain  the  inviscid  properties  along  the  boundary- layer 
edge  away  from  the  stagnation  point.  The  ideal  gas  (7  =  1.4) 
blunt  body  and  characteristics  solutions  provided  the  surface 
pressure  distribution  used  in  the  streamtube  expansions.  It 
should  be  noted  that  the  choice  of  inviscid  nonequilibrium 
stagnation  conditions  rather  than  equilibrium  stagnation  con¬ 
ditions  was  made  on  practical  rather  than  theoretical  grounds. 
In  most  previous  studies,  it  has  been  assumed  that  the  cor¬ 
rect  inviscid  stagnation  conditions  were  thermodynamic  equi¬ 
librium.  This  model  is  theoretically  correct  if  the  boundary 
layer  is  sufficiently  thin  with  respect  to  the  shock  layer 
thickness  in  the  stagnation  region.  However,  if  the  boundary 
layer  is  thick  with  respect  to  the  shock  layer,  the  condi¬ 
tions  at  the  "edge”  will  not  be,  in  general,  the  equilibrium 
stagnation  conditions.  In  attempting  to  account  for  this 
departure  from  chemical  equilibrium  in  the  stagnation  region, 
the  inviscid  nonequilibrium  blunt  body  solution  method  of 
Curtis  and  Strom  was  used  to  obtain  the  inviscid  stagnation 
conditions.  The  stagnation  streamline  is  not  included  in  the 
blunt  body  solution  method  of  Ref.  12,  and  the  stagnation 
point  data  were  obtained  from  extrapolations  of  the  solution 
from  sQ  =  0.1  to  s  =  0.  The  nonequilibrium  streamtube  expan¬ 
sion  was  begun  at  s  =  sQ  with  the  initial  data  from  the  blunt 
body  solution  at  the  same  location.  Attempts  to  obtain  non¬ 
equilibrium  solutions  over  the  conical  afterbody  were  not 
successful,  and  therefore  the  pressure  distributions  used  in 
the  streamtube  expansions  were  obtained  from  the  ideal  gas 
(7  =  1.4)  blunt  body  and  characteristics  solutions.  It 
should  be  noted  that  the  nonequilibrium  blunt  body  solution 
was  completely  inviscid  and  did  not  include  the  boundary- layer 
option  described  by  Curtis  and  Strom  (Ref.  12). 

The  effects  of  chemical  nonequilibrium  on  the  inviscid 
stagnation  conditions  were  considered.  Stagnation  point 
boundary- layer  calculations  were  made  using  both  chemical 
equilibrium  and  nonequilibrium  inviscid  conditions.  The  ef¬ 
fects  on  the  stagnation  heat- transfer  rate  were  considered, 
and  the  results  are  presented  below. 


2.2  VISCOUS  FLOW  FIELD  MODELS 

Most  of  the  calculations  presented  in  this  report  were 
made  using  a  binary  perfect  gas,  laminar  boundary- layer  model 
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including  transverse  curvature  and  displacement  effects. 
Limited  calculations  are  presented  from  an  ionizing  air,  thin 
laminar  boundary- layer  theory  and  a  thin  viscous  shock-layer 
model.  The  theoretical  models  are  described  below,  and  the 
conditions  used  in  the  calculations  are  given  in  a  later  sec¬ 
tion. 

2.2.1  Binary  Perfect  Gas  Boundary-Layer  Model 

The  method  of  Jaffe,  Lind,  and  Smith  (Ref.  5)  was  used 
to  solve  the  laminar  boundary-layer  equations  including 
transverse  curvature  effects  for  a  binary  mixture  of  helium 
or  argon  and  air.  All  gases  were  treated  as  perfect  and, 
for  these  calculations,  the  thermodynamic  and  transport  pro¬ 
perties  were  extended  as  described  in  Appendix  III. 

It  was  not  possible  to  match  all  of  the  wind  tunnel 
free-stream  and  stagnation  conditions  with  the  perfect  gas 
data  in  the  binary  gas  boundary- layer  program.  The  procedure 
used  was  to  compute  the  stagnation  enthalpy  from  H*,  =  cp  T*, 

+  1/2  m2  -yRToo  with  cp  =  6011.6  ft2/sec2-°R,  Tro  -  274°R,  = 

13.2,  7  =  1.4,  and  R  =  1717.6  ft2/sec2-°R.  The  inviscid 
normal  shock  stagnation  temperature  Tq  was  then  computed  from 
the  thermodynamic  properties  in  the  binary  gas  boundary- layer 
program,  and  then  an  effective  specific  heat  ratio  7  was 
computed  from  7/(7  -  1)  =  H^/CRT').  The  ideal  gas  (7  -=  7) 
normal  shock  stagnation  pressure  was  calculated  from  the 
Rayleigh  pitot  formula.  The  results  of  this  procedure  are 
given  below: 


Ames  Tunnel  Conditions 

Binary  Gas 

HmJ  ft2/sec2 

6.509  x  107 

5.9047  x  10 

Mco 

00 

13.2 

13.2 

U  ,  ft/sec 

00 

10,850 

10,680 

T  ,  °R 

274 

290 

P00,  lbf/ft2 

0.1944 

0.1944 

7  >  7 

1.4 

1.313 

T',  °R 

— 

8204 

p',  lbf/ft2 

40.4 

41.65 

4 


AEDC-TR-68-214 


The  free-stream  Mach  number  and  pressure  were  matched  and 
the  remaining  quantities  were  calculated  from  ideal  gas  re¬ 
lations.  The  input  conditions  at  the  edge  of  the  boundary 
layer  were  obtained  from  an  ideal  gas  (7  =  7)  one- dimensional 
expansion  using  the  ideal  gas  (7  =  1.4)  blunt  body  and  charac¬ 
teristics  solutions.  The  effects  of  the  mismatch  in  experi¬ 
mental  and  theoretical  conditions  due  to  the  ideal  gas  model 
chosen  are  difficult  to  determine;  however,  the  results  of 
previous  studies  such  as  Whitfield  and  Griffith  (Ref.  1)  and 
Griffith  and  Lewis  (Ref.  14)  where  both  perfect  gas  and 
equilibrium  air  calculations  were  made  indicated  small  ef¬ 
fects  due  to  variations  in  boundary- layer  edge  conditions  with 
gas  properties. 

The  procedure  used  to  calculate  the  effects  of  boundary- 
layer  displacement  is  called  a  first-order  treatment.  The 
inviscid  outer  flow  is  perturbed  by  the  boundary-layer  dis¬ 
placement  thickness  including  the  effects  of  mass  transfer 

6*  through 
m 

5*  e 

m  6  x 

/  Peuer  dy  =  /  (peue  -  Pu>r  dy  +  /  Pwvwrw  dx  U) 

o  o  o 

where,  without  mass  transfer,  the  displacement  thickness  6* 

-  <5*,  which  is  the  axisymmetric  form  of  the  classical  two- 
dimensional  displacement  thickness.  It  should  be  noted  that 
this  is  not  the  definition  proposed  by  Fannelop  (Ref.  15) 
since  he  neglected  higher-order  terms  which  are  significant 
under  low  density  conditions. 

The  inviscid  outer  flow  sees  an  effective  body  defined 
by 

r  =  r  +  6*  cos  a  (2) 

eff  w  m 

The  effect  of  displacement  was  obtained  by  iterating  the 
inviscid  outer  flow  with  the  viscous  inner  flow  until  there 
was  negligible  change  in  the  effective  body  pressure  distri¬ 
bution.  Consistent  with  classical  boundary- layer  theory,  the 
normal  pressure  gradient  was  ignored,  and  the  effects  of 
shock-generated  external  vorticity,  slip,  and  temperature 
jump  were  also  not  considered  in  the  present  study.  It  has 
been  shown  by  Lewis  (Ref.  6)  and  Marchand,  Lewis,  and  Davis 
(Ref.  11)  that  strong  coupling  effects  can  exist  between 
higher-order  boundary-layer  effects  which  can  significantly 
affect  wall  shear  stress  and  heat- transfer  rates.  However, 
since  numerical  methods  have  not  been  developed  to  consider 
binary  gas  or  ionizing  air  boundary  layers  including  all  the 
higher-order  effects,  these  first-order  results  including 


5 


AEDC-TR-68-214 


tranverse  curvature  and  displacement  should  provide  useful 
estimates  for  engineering  purposes. 

The  boundary  conditions  on  the  species  continuity  equa¬ 
tions  originally  proposed  by  Jaffe,  Lind,  and  Smith  (Ref.  5) 
were  later  changed  (Ref.  16)  to  agree  with  the  proposal  of 
Baron  (Ref.  17).  It  was  originally  proposed  in  Ref.  5  that 
the  concentration  of  the  injected  gas  be  100  percent  at  the 
wall  and,  of  course,  zero  at  the  outer  edge.  Baron  proposed 
that  the  correct  wall  boundary  condition  be  zero  net  flux  of 
the  outer  flow  component  (air).  Using  Baron's  idea,  the  wall 
boundary  conditions  on  both  the  specie  conservation  equation 
(Zw)  and  the  energy  equation  (gw)  are  unknown  a  priori.  One 
must  therefore  solve  for  the  unknown  boundary  conditions  by 
iteration  using  the  condition  of  zero  net  flux  of  air  at  the 
wall  to  determine  the  injected  gas  specie  mass  fraction  and 
the  mixture  enthalpy  at  the  wall.  A  numerical  scheme  was 
developed  by  Mayne,  Gilley,  and  Lewis  (Ref.  4)  based  on 
Baron's  proposal,  and  that  method  was  used  in  the  present 
calculations . 

2.2.2  Viscous  Nonequilibrium  Flow  Field  Analysis 


The  determination  of  the  viscous  flow  field  near  a  body 
surface  requires  a  detailed  calculation  of  the  boundary- layer 
equations.  This  task  was  accomplished  through  the  numerical 
solution  of  these  governing  relations  which  provide  convec¬ 
tive  heat  flux  results  as  well  as  velocity,  temperature, 
density,  and  specie  concentration  profiles  across  the  boundary 
layer.  Furthermore,  the  effects  of  inert  gas  injection  into 
the  boundary  layer  have  been  included  so  as  to  permit  study  of 
the  transpiration  cooling  problem  under  nonequilibrium  condi¬ 
tions  . 

The  computer  program  used  to  solve  the  boundary- layer 
equations  employs  a  highly  modified  version  of  the  method 
developed  by  Blottner  (Ref.  8)  to  study  nonequilibrium,  non¬ 
similar,  laminar  boundary- layer  flows  over  both  blunt  and 
sharp-nosed  bodies  of  revolution  as  well  as  sharp  plates  and 
wedges.  All  basic  portions  of  Blottner 's  analysis  have  been 
retained,  e.g.,  consideration  of  the  diffusive  mass  flux  for 
a  multicomponent  mixture  of  perfect  gases  and  expressing  the 
net  mass  rate  of  production  for  each  specie  in  a  form  espe¬ 
cially  suited  for  numerical  computation.  For  further  details 
on  the  theoretical  formulation  of  the  governing  equations, 
see  Blottner  (Refs.  18  and  19)  and  Lenard  (Ref.  20). 

In  order  to  examine  the  effects  of  inert  gas  injection 
into  a  nonequilibrium  boundary  layer,  the  wall  boundary  con¬ 
ditions  for  the  governing  specie  conservation  equations  must 
be  modified  to  allow  mass  transfer  through  a  porous  surface. 
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Following  Chung  (Ref.  21)  and  Moore  and  Lee  (Ref.  22),  the 
solid  boundary  (wall)  was  viewed  as  a  rigid  semipermeable 
membrane  which  does  not  permit  the  air  species  to  flow  into 
the  pores  but  does  allow  the  inert  species  to  flow  into  the 
boundary  layer.  Although  chemically  inert,  many  surfaces 
can  act  as  a  catalyst  for  the  recombination  of  oxygen  and 
nitrogen  atoms  in  a  dissociated  mixture  so  that  a  catalytic 
wall  boundary  condition  must  also  be  considered.  Moore 
(Ref.  23),  in  an  appendix  to  his  thesis,  discusses  and 
formulates  these  boundary  conditions  for  two  different  cases 
with  inert  gas  injection,  one  in  which  the  wall  is  taken  to 
be  completely  noncatalytic  with  respect  to  atom  recombination 
and  the  other  in  which  local  chemical  equilibrium  is  assumed 
to  exist  at  the  gas-solid  interface.  In  the  absence  of  mass 
transfer  these  two  boundary  conditions  reduce  to  the  non¬ 
catalytic  and  equilibrium  wall  conditions,  respectively, 
formulated  by  Blottner  and  Lenard  (Ref.  24).  Hence,  Moore's 
formalism  of  the  wall  boundary  conditions  for  the  species  may 
be  considered  general  and  used  both  with  and  without  mass 
transfer  effects;  these  boundary  conditions  were  used  in  the 
present  study  by  adding  them  to  the  computer  program.  The 
outer  edge  boundary  conditions  were  obtained  from  the  non¬ 
equilibrium  blunt  body  and  streamtube  expansion  as  discussed 
previously.  In  addition,  the  usual  conditions  of  no  slip 
and  temperature  jump  at  the  surface  and  a  prescribed  constant 
wall  temperature  (540°R)  were  made. 

In  addition  to  the  wall  and  outer  edge  boundary  condi¬ 
tions,  initial  conditions  are  required  to  completely  determine 
the  resulting  downstream  solution.  In  order  to  obtain  the 
correct  initial  conditions,  the  system  of  boundary-layer 
equations  (which  are  nonlinear  partial  differential  equations) 
was  reduced  to  the  limiting  form  of  ordinary  nonlinear  differ¬ 
ential  equations  with  two-point  boundary  conditions  at  the 
stagnation  point.  Such  two-point  boundary  value  problems  are 
extremely  difficult  to  solve  using  the  more  common  "shooting" 
method  since  so  many  specie  equations  are  involved.  Hence, 
a  finite-difference  method  was  developed  such  that  the  two- 
end  conditions  were  always  satisfied;  the  fundamentals  of 
this  technique  may  be  found  in  Chapter  VII  of  the  book  on 
numerical  analysis  by  Conte  (Ref.  25).  More  details  of  this 
finite-difference  technique  applied  to  the  solution  of  non¬ 
linear  equations  may  be  found  in  the  papers  by  Holt  (Ref.  26) 
and  Fay  and  Kaye  (Ref.  27).  Basically,  the  nonlinear  ordinary 
differential  equations  were  written  in  linearized  finite- 
difference  form  so  that  the  resulting  linear  algebraic  equa¬ 
tions  are  of  tridiagonal  form  and  a  sweeping  method  of  solu¬ 
tion  was  used.  Iteration  was  then  applied  until  successive 
solutions  (f'(Ti),  0(tj),  and  C^(t|))  were  sufficiently  close 
to  each  other.  This  procedure  was  found  to  converge  rapidly 
for  all  cases  in  the  present  study. 
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Using  the  above-obtained  initial  conditions,  solutions 
to  the  complete  nonsimilar  boundary- layer  equations  were  ob¬ 
tained  in  a  step-by-step  procedure  marching  downstream  using 
the  implicit  finite-difference  method  by  Blottner  which  was 
mentioned  previously.  No  stability  problems  were  encountered, 
either  with  or  without  mass  transfer,  for  the  conditions 
studied. 

Considerations  have  also  been  given  to  the  possibility 
that  a  boundary-layer-type  analysis  may  not  be  valid  under 
the  conditions  of  the  present  study.  It  may  be  that  non¬ 
equilibrium  and  viscous  effects  have  equal  importance  over 
the  entire  shock  layer,  especially  in  the  nose  region  of  the 
body.  In  order  to  examine  this  possibility,  an  analysis  has 
been  performed  where  the  entire  shock  layer  at  the  stagnation 
point  was  treated  as  viscous.  As  shown  by  Cheng  (Ref.  28), 
Chung  (Ref.  29),  Inger  (Ref.  30),  Howe  and  Viegas  (Ref.  31), 
and  Lee  and  Zierten  (Ref.  32),  the  governing  equations  for  such 
an  analysis  are  identical  in  form  with  the  boundary- layer 
equations  under  the  restrictions  that  the  shock  layer  be  thin, 
i.e.,  A  «  1,  with  all  curvature  effects  neglected.  Recent 
work  by  Howe  and  Sheaf fer  (Ref.  33),  Goldberg  (Ref.  34), 
Goldberg  and  Scala  (Ref.  35),  and  Chen,  et  al .  (Ref.  36),  has 
further  shown  the  validity  of  such  an  analysis  for  cases  with 
mass  injection  at  the  wall. 

The  basic  thin,  hypersonic  shock  layer  model  as  formu¬ 
lated  by  Cheng  assumes  a  strong  bow  shock  wave  of  negligible 
thickness  that  is  concentric  with  the  body  at  the  stagnation 
line.  Across  the  shock,  the  static  pressure  and  normal 
velocity  components  are  discontinuous  according  to  the  usual 
Rankine-Hugoniot  relations;  in  addition,  the  so-called 
"shock  slip"  effects  (discontinuities  in  composition,  tan¬ 
gential  velocity,  and  total  enthalpy)  are  allowed  immediately 
behind  the  shock.  Between  the  shock  and  the  body  there  is 
assumed  a  thin  (A  «  1)  layer  of  continuum,  viscous  reacting 
gas  flow.  In  this  layer,  the  influence  of  body  curvature 
and  surface  slip  phenomena  are  also  neglected  in  comparison 
to  the  shock  layer  vorticity  and  "shock  slip"  effects,  pro¬ 
vided  the  body  surface  is  highly  cooled.  An  important  con¬ 
sequence  of  this  thin,  hypersonic  shock  layer  model  is  that 
the  governing  Navier-Stokes  equations  reduce  to  the  familiar 
parabolic  partial  differential  equations  of  boundary- layer 
theory  as  mentioned  previously. 

For  application  in  the  present  study,  a  modification  of 
the  above  model  was  adopted.  Across  the  shock,  the  static 
pressure,  enthalpy,  and  normal  velocity  component  were  taken 
to  be  discontinuous  according  to  the  Rankine-Hugoniot  condi¬ 
tions  for  a  real  gas.  However,  through  the  shock  the  chem¬ 
ical  composition  was  taken  as  frozen  at  the  free-stream 
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conditions.  Furthermore,  stagnation  or  total  enthalpy  was 
taken  as  constant  across  the  shock.  Hence,  none  of  the  so- 
called  "shock  slip"  effects  was  considered  in  the  current 
work.  Other  than  this,  all  other  features  of  the  thin, 
hypersonic  shock  layer  were  retained  including  the  centrif¬ 
ugal  effects  due  to  a  normal  pressure  gradient  across  the 
layer.  Between  the  shock  and  the  body  the  flow  was  con¬ 
sidered  to  be  in  chemical  nonequilibrium  with  the  shock 
standoff  distance  determined  as  part  of  the  solution.  Ef¬ 
fects  of  inert  gas  injection  into  the  shock  layer  were  also 
included  for  both  noncatalytic  and  catalytic  wall  conditions. 
The  program  itself  is  basically  the  same  as  the  stagnation 
point  boundary- layer  program  described  previously  since  the 
only  modifications  were  to  permit  determination  of  the  shock 
standoff  distance  and  the  inclusion  of  the  normal  pressure 
gradient  terms.  The  manner  of  iterative  solution  was  exactly 
the  same  as  in  the  boundary- layer  case.  Future  plans  include 
the  addition  of  "shock  slip"  effects  to  this  program  since 
only  with  this  capability  can  the  low  Reynolds  number  regime 
of  high  altitude  hypersonic  flight  be  explored. 

SECTION  III 

CONDITIONS  USED  IN  THE  CALCULATIONS 


The  free-stream  and  stagnation  conditions  used  in  the 
binary  gas  boundary- layer  calculations  are  given  in  Table  I 
(Appendix  II) .  The  injection  rates  were  assumed  constant 
and  are  given  below: 


Injected  Gas 

Argon  (pwvw,  lbf-sec/ft3) 
Helium  (pwvw,  lbf-sec/ft3) 

Argon  (pwvw/poouoo) 

Helium  (p^Vp^uJ 


Binary  Gas 

11.82  x  10"5 
3.08  x  10~5 

0.0284 

0.0074 


Ionizing  Air 

11.82  x  10-5 
3.08  x  10"5 

0.0268 

0.0070 


It  will  be  shown  later  that  the  higher  injection  rate  of  the 
argon  had  a  significant  effect  on  the  numerical  results. 

The  porous-wall  temperature  was  assumed  constant  at 
Tw  =  540°R . 
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The  wall  and  outer-edge  specie  mass  fractions  for  the 
ionizing  air  boundary  layer  and  shock  layer  are  given  in 
Tables  II  through  IV.  Where  computed,  both  equilibrium 
catalytic  wall  and  noncatalytic  wall  conditions  are  given. 

The  equilibrium  wall  conditions  were  computed  from  the 
equilibrium  composition  at  the  wall  (or  outer  edge)  pressure 
and  wall  temperature.  With  injection  the  added  constraint 
of  zero  net  flux  of  the  air  species  at  the  wall  was  imposed. 
The  frozen  normal  shock  specie  mass  fractions  for  the  thin 
viscous  shock  layer  are  equivalent  to  the  free-stream  mass 
fractions . 

For  comparison  of  the  effects  of  equilibrium  and  non¬ 
equilibrium  inviscid  outer  edge  conditions,  no- injection  thin 
boundary- layer  calculations  were  made  with  both  sets  of  con¬ 
ditions  and  the  inviscid  outer  edge  mass  fractions  are  given 
in  Table  II. 

The  chemical  reaction  rate  constants  used  in  the  non¬ 
equilibrium  boundary- layer  and  shock-layer  calculations  are 
given  in  Table  V. 


SECTION  IV 

RESULTS  AND  DISCUSSION 


4.1  BINARY  GAS  RESULTS 

The  results  from  the  binary  perfect  gas,  laminar  boundary- 
layer  calculations  including  the  effects  of  displacement  and 
transverse  curvature  are  shown  in  Figs.  2-7  for  the  no¬ 
injection  conditions,  in  Figs.  8-14  for  the  argon- injection 
conditions,  and  in  Figs.  15-21  from  the  helium-injection 
calculations . 


The  effects  of  displacement  on  all  quantities  plotted  can 
be  seen  by  comparison  of  the  last  iteration  (e.g.  second 
iteration  without  injection)  value  with  the  zeroth- iteration 
results.  The  convergence  of  the  solutions  with  displacement 
can  be  seen  by  comparison  of  (n  +  1)-  with  n-iteration  re¬ 
sults.  From  the  results  shown  in  Figs.  8-13,  it  is  clear 
that  the  convergence  of  the  calculations  for  the  displacement 
effect  with  argon  injection  were  much  slower  (i.e.,  more 
iterations  were  required  for  negligible  change  between  iter¬ 
ations)  than  either  the  helium*  or  no-injection  calculations. 


Because  of  an  error  pointed  out  to  the  authors  by  Dr. 

N.  A.  Jaffe,  the  helium- inject ion  calculations  were  corrected 
and  only  the  zeroth  and  final  interation  results  are  shown. 

The  error  had  negligible  effect  on  the  argon-injection  results. 
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The  higher  injection  rate  used  in  the  argon  calculation 
caused  the  slower  convergence.  It  should  be  noted  that  dif¬ 
ficulty  was  encountered  in  making  the  first-iteration  argon 
injection  calculation  at  s  =  4.32  which  caused  the  calcula¬ 
tion  to  be  terminated  at  that  point.  The  effective  body  data 
were  extrapolated  over  the  remainder  of  the  body,  and  the 
calculation  was  continued  for  the  second  iteration.  This 
difficulty  with  high  injection  rates  was  encountered  in  the 
earlier  sharp  cone  studies  of  Lewis,  Marchand,  and  Little 
(Ref.  2)  and  Mayne,  Gilley,  and  Lewis  (Ref.  4),  and  the 
extrapolation  procedure  used  in  the  present  study  was  also 
used  successfully  in  the  previous  studies. 

The  primary  interest  in  the  present  investigation  was  in 
the  predictions  of  heat- transfer  and  pressure  distributions 
over  the  cone  with  and  without  injection.  The  he at- transfer 
distributions  are  shown  in  dimensional  (Btu/ft2-sec)  and  two- 
dimensionless  forms  (4/p^U^  and  q/q(0)).  The  data  are  not 
shown  in  the  form  of  the  Stanton  number  since  the  enthalpy 
of  the  gas  mixture  at  the  wall  is  a  function  of  the  concen¬ 
tration  of  the  injected  gas  at  the  wall  which  is  unknown 
a  priori.  Also  the  driving  potential  for  heat  transfer  is 
not  simply  the  enthalpy  difference  when  injection  or  chemical 
reactions  are  considered  but  is  a  more  complicated  function. 
An  alternate  and  somewhat  preferable  nondime ns ional  heat- 
transfer  coefficient  is  q/p  U2  which  is  useful  in  scaling 
results  for  small  variations  in  free-stream  conditions.  The 
commonly  used  ratio  q/q(0)  is  also  given. 

The  mass  fractions  of  the  injected  gases  (argon  and 
helium)  at  the  wall  are  shown  in  Figs.  14  and  21.  The  high 
uniform  injection  rate  of  the  argon  increased  the  mass  frac¬ 
tion  of  argon  from  0.14  at  the  stagnation  point  to  0.85  at 
the  end  of  the  body,  whereas  the  variation  in  helium  mass 
fraction  increased  from  0,025  at  s  =  0  to  0.44  at  s  =  5.10. 

It  has  been  observed  previously  that  blowoff  of  the  boundary 
layer  (i.e.,  where  tw  ~  f ^  0)  is  related  to  the  concentra¬ 

tion  of  the  injected  gas  approaching  100  percent  when  uniform 
surface  mass  transfer  is  considered. 

The  effects  of  injection  on  the  heat-transfer  are  shown 
in  Fig.  22  with  a  plot  of  the  ratio  of  heat  transfer  with 
helium  or  argon  injection  to  the  no-injection  heat-transfer 
rate.  Both  displacement  and  transverse  curvature  effects 
were  included  in  the  injection  and  no-injection  calculations 
with  the  final  iteration  results  used  in  all  cases.  The 
higher  injection  rate  of  the  argon  is  seen  to  substantially 
reduce  the  heat- transfer  rate  in  comparison  with  the  helium 
injection  results. 


11 


AEDC-TR-68-214 


Attention  is  directed  to  the  increase  in  heat- transfer 
in  the  stagnation  region  (s  <  0.85)  with  helium  injection 
over  the  no-injection  heat-transfer  rate  as  shown  in  Fig.  22. 
Some  of  the  quantities  affecting  the  wall  heat- transfer  rate 


are  given  below: 

Helium  Injection 

No  Injection 

V  °R 

540 

540 

pw,  slugs/ft3 

3.93216  x  10“5 

4.49397  x  10“5 

p.  ,  lbf-sec/ft2 

W 

3.5938  x  10“7 

3.50729  x  10"7 

2  2 

c_  ,  ft  /sec  -R 

Pw 

6593 

6019 

TY 

Pr 

w 

0.60447 

0.73725 

Le 

w 

2.8649 

1.00 

kw,  lbf/sec-°R 

3.9198  x  10"3 

2.8634  x  10"3 

Zw 

0.97705 

1.00 

V 

0.003347 

0.0 

gw 

0.060274 

0.055026 

gw 

0.205042 

0.2158078 

q,  f t-lbf/f t2-sec 

38,146 

36,205 

(dT/dy)w 

0.81  x  106 

1.06  x  106 

(dZ/dy)w 

1.466 

0.0 

The  heat-transfer  rate  to  the  wall  with  or  without  mass 
transfer  is  shown  in  Appendix  IV  to  be  =  (k  dT/dy)w  for  a 
nonreacting  binary  gas  mixture.  For  this  condition,  mass 
transfer  can  only  affect  the  wall  heat  transfer  through  the 
mixture  thermal  conductivity  and  temperature  gradient.  Con¬ 
sidering  the  ratio 

^He  inj  =  (k  5T/Sy)He  inj,w 

^no  inj  ^k  ^T/dy)no  inj,w 

and  the  data  given  above,  one  finds  the  ratio  to  be  1.05  in 
good  agreement  with  the  results  shown  in  Fig.  22.  For  the 
case  of  helium  injection,  it  is  thus  seen  that  the  increase 
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in  thermal  conductivity  of  the  gas  mixture  at  the  wall  is 
greater  than  the  reduction  in  the  temperature  gradient  at  the 
wall.  This  interesting  result  was  also  found  by  Anfimov 
(Ref.  37)  for  injection  of  low  molecular  weight  gases  (He  and 
H2)  into  a  nonreacting  binary  gas  mixture.  In  agreement  with 
the  results  given  later  in  this  report,  he  did  not  find  the 
increase  in  stagnation  point  heat  transfer  when  He  was  in¬ 
jected  into  a  reacting  air  mixture.  This  phenomenon  appears 
to  be  confined  to  low  injection  rates  of  low  molecular  weight 
gases  into  nonreacting  binary  gas  mixtures.  It  will  be  of 
special  interest  to  compare  these  results  with  experimental 
heat-transfer  data  to  determine  the  magnitude  and  direction 
of  the  change  in  surface  heat-transfer  rate  with  low  injec¬ 
tion  rates  of  helium. 


4.2  NONEQUILIBRIUM  FLOW  FIELD  RESULTS 
4.2.1  Stagnation  Point 


a 

It  is  instructive  first  to  examine  the  stagnation  point 
solutions.  Information  concerning  the  no  mass  injection  con¬ 
dition  is  shown  in  Table  II.  These  basic  sets  of  solutions 
are  shown:  a  boundary- layer  model  which  used  the  nonequilib¬ 
rium  inviscid  outer  edge  specie  concentrations,  a  boundary- 
layer  model  which  used  equilibrium  outer  edge  specie  concen¬ 
trations  as  computed  from  equilibrium  considerations  at  the 
given  pQ  and  Tq  (note  that  these  are  perfect  gas  values)  with 

reference  mass  fractions  *=  0.23456  and  C,lre:f^  =  0.76544, 

O2  N2 

and  a  thin  viscous  shock-layer  model  using  the  conditions 
behind  the  chemically  frozen  (free-stream  composition)  normal 
shock.  Note  the  large  difference  in  outer  edge  specie  con¬ 
centrations  between  the  boundary- layer  nonequilibrium  and 
equilibrium  cases;  the  equilibrium  values  for  oxygen  are 
highly  dissociated  and  slightly  dissociated  for  nitrogen, 
whereas  the  nonequilibrium  values  are  essentially  frozen  at 
their  free-stream  level  for  which  the  oxygen  was  slightly 
dissociated.  This  rather  striking  difference  becomes  strongly 
reflected  in  the  stagnation  point  heat-transfer  results  for 
the  fully  catalytic  wall  condition.  However,  before  examin¬ 
ing  this,  a  review  of  various  wall  conditions  and  their 
interpretation  is  in  order. 

As  used  here,  the  fully  catalytic  wall  (FCW)  condition 
is  taken  to  mean  that  the  wall  is  perfectly  catalytic  with 
respect  to  atom  recombination  so  that  only  the  basic  molecular 
species  oxygen  and  nitrogen  exist  at  the  surface.  Further, 
it  is  assumed  that  the  values  for  oxygen  and  nitrogen  com¬ 
position  at  the  wall  are  exactly  the  same  as  the  outer  edge 
reference  values  used  in  the  equilibrium  composition  calculation. 
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Such  is  purely  an  assumption  since  in  the  case  of  a  multi- 
component  gas  the  basic  element  mass  fractions  may  change  due 
to  multicomponent  diffusion  across  a  nonisothermal  layer. 

This  assumption  can  only  be  justified  a  posteriori  by  examin¬ 
ing  the  complete  multicomponent  formalism  for  equilibrium 
catalytic  wall  (ECW)  boundary  conditions.  Such  has  been  done 
for  the  thin  viscous  shock  layer  case,  and  the  results  are 
shown  in  Table  VI.  It  is  seen  that  the  fully  catalytic  wall 
assumption  is  certainly  valid  for  this  case;  repetition  of 
the  equilibrium  wall  calculations  for  the  boundary-layer 
cases  established  the  same  conclusion  with  respect  to  the 
fully  catalytic  wall  assumption  for  these  conditions.  It 
should  be  pointed  out  that  the  validity  of  this  assumption 
is  due  mainly  to  the  extremely  cold  wall  temperature  (540°R) 
used  throughout  this  analysis.  x  As  the  wall  temperature  is 
increased,  the  fully  catalytic  wall  assumption  will  lose  its 
validity  and  recourse  must  be  taken  to  the  multicomponent 
equilibrium  wall  formalism. 

An  alternate  boundary  condition  for  the  species  is  the 
noncatalytic  wall  (NCW)  restriction.  Here  no  atom  recombina¬ 
tion  is  allowed  on  the  surface,  and  thus  the  atoms  tend  to 
"pile  up"  near  the  wall,  creating  a  "blanketing"  effect.  The 
model  surface  will  exhibit  characteristics  somewhere  between 
the  two  limiting  extremes  described  above,  i.e.,  a  finite 
surface  catalyticity  with  respect  to  atom  recombination. 

This  condition,  which  is  difficult  to  treat  due  to  lack  of 
knowledge  concerning  surface  catalyticity  effects,  has  not 
been  considered  in  the  present  analysis  since  the  fully 
catalytic  (or  equilibrium)  wall  and  the  noncatalytic  wall 
form  an  upper  and  lower  bound,  respectively,  on  the  heat- 
transfer  rate  with  respect  to  atom  recombination  on  the  surface. 

Since  the  noncatalytic  wall  condition  allows  no  atom  re¬ 
combination  on  the  surface,  the  wall  heat  transfer  can  be 
materially  reduced  for  cases  involving  high  dissociation  (and 
ionization)  levels  in  the  gas.  An  excellent  example  of  this 
effect  may  be  seen  in  the  boundary-layer  solution  with  equi¬ 
librium  outer  edge  conditions.  Note  the  factor  of  two 
difference  in  the  stagnation  point  heat  flux;  this  is  due  to 
the  recombination  of  atoms  on  the  surface  in  the  fully 
catalytic  wall  case  as  opposed  to  no  recombination  on  the 
surface  in  the  noncatalytic  wall  case  under  conditions  of  high 
dissociation  levels  at  the  outer  edge  of  the  boundary  layer. 

Such  does  not  occur  in  the  boundary-layer  case  with  a  non¬ 
equilibrium  outer  edge  because  here  the  dissociation  level  is 
very  low,  which  is  reflected  in  a  very  modest  increase  in 
stagnation  point  heat  flux  for  the  fully  catalytic  wall  con¬ 
dition  in  comparison  with  the  noncatalytic  wall  conditions. 

It  should  be  further  noted  that  the  flow  is  essentially 
frozen  between  the  outer  edge  and  the  wall  for  the  non¬ 
catalytic  wall  condition,  and  this  will  be  shown  later  in  the 
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discussion.  The  important  fact  to  be  gained  from  the  above 
is  that  the  assumption  on  the  outer  edge  condition  (either 
equilibrium  or  nonequilibrium)  can  be  controlling  with  respect 
to  the  resulting  stagnation  point  heat  flux.  Hence,  one 
should  try  to  determine  the  best  choice  for  the  outer  edge 
condition  from  some  independent  source  if  possible.  The  thin 
viscous  shock  layer  analysis  provides  such  a  tool. 

Since  the  thin  viscous  shock  layer  (TVSL)  analysis  treats 
the  flow  field  as  fully  viscous  and  in  chemical  nonequilibrium 
between  the  shock  and  the  body,  there  is  no  arbitrary  choice 
as  to  the  location  of  some  "edge"  (with  its  corresponding 
conditions)  since  the  shock  standoff  is  now  obtained  as  part 
of  the  solution.  The  only  restriction  made  is  that  the  normal 
shock  crossing  be  chemically  frozen  with  translational,  rota¬ 
tional,  and  vibrational  modes  in  equilibrium  immediately  be¬ 
hind  the  shock.  Reference  to  Table  II  shows  that  the  stagna¬ 
tion  point  heat  fluxes  are  in  good  agreement  with  the  non¬ 
equilibrium  outer  edge  boundary- layer  results  for  both  wall 
conditions.  This  agreement  may  be  interpreted  as  evidence 
that  perhaps  the  outer  edge  of  the  boundary- layer  analysis 
should  properly  be  treated  as  a  nonequilibrium  condition.  It 
should  be  noted,  however,  that  vorticity  is  included  in  the 
viscous  shock  layer  model  (increasing  the  heat-transfer  rate) 
and  not  in  the  classical  boundary- layer  model.  Thus,  the 
good  agreement  in  wall  heat-transfer  between  the  two  models 
must  be  viewed  with  caution.  Of  course,  one  may  even  question 
the  applicability  of  a  thin  viscous  shock  layer  treatment  for 
the  present  set  of  conditions.  In  order  to  verify  this 
analysis,  the  shock  Reynolds  number,  Res,  must  be  examined. 

For  the  particular  set  of  free-stream  conditions  shown  in 
Table  I,  Res  =  149.93  based  on  body  nose  radius;  the  corre¬ 
sponding  density  ratio,  e,  across  the  frozen  normal  shock 
was  e  =  0.1389  with  static  pressure  ps  =  41.471  lbf/ft2  and 
static  temperature  Ts  =  8069. 2°R  immediately  behind  the  shock. 
As  shown  by  Cheng  (Ref.  28),  the  conditions  for  validity  of 
a  thin  viscous  shock  layer  model  are  (l/€Res)  «  1,  A  =  0(e) 

«  1  and  eM£  »  1 •  For  the  present  conditions 

l/eRe_  =  1/(149.93  x  0.1389)  =  0.0480  «  1 

A  =  0(e)  =  0(0.1389)  «  1 

eM^  =  (0.1389) (13. 2)2  =  24.2019  »  1 

which  satisfy  these  criteria  and,  hence,  justify  a  thin  vis¬ 
cous  shock  layer  analysis.  Thus,  it  is  believed  that  the 
thin  viscous  shock  layer  is  an  applicable  model  upon  which 
to  assess  the  validity  of  a  boundary- layer- type  model  by  com¬ 
parison  of  solutions'. 
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From  Table  II  it  is  seen  that  the  inviscid  shock  stand¬ 
off  A  is  0.1091  for  the  inviscid  nonequilibrium  solution, 

0.1095  for  the  inviscid  perfect  gas  solution,  and  0.0753  or 
0.0758  for  the  thin  viscous  shock  layer  fully  catalytic  and 
noncatalytic  wall,  respectively.  The  reduction  in  shock 
layer  thickness  for  the  thin  viscous  shock  layer  cases  is  due 
to  the  large  density  gradients  near  the  cold  wall  which  tend 
to  thin  the  layer.  For  the  boundary  layer  with  nonequilibrium 
outer  edge,  the  boundary- layer  thickness  6  (defined  to  be  the 
y  location  where  f '  =  0.995)  is  seen  to  be  0.0695  and  0.0696 
for  the  fully  catalytic  and  noncatalytic  wall,  respectively. 
Hence,  the  viscous  boundary- layer  occupies  about  65  percent 
of  the  inviscid  shock  layer;  note  further  that  the  boundary- 
layer  thickness  is  slightly  less  than  the  thin  viscous  shock- 
layer  thickness.  All  of  these  factors  contribute  to  the 
possibility  that  a  boundary-layer-type  analysis  is  indeed  a 
valid  model  for  the  present  body  and  set  of  conditions. 

With  respect  to  mass  addition  at  the  stagnation  point, 
Tables  III  and  IV  show  results  for  the  injection  of  argon  and 
helium,  respectively.  First  note  that  the  injected  mass  flux 
for  helium  is  roughly  one  order  or  magnitude  less  than  for 
argon;  such  has  an  important  effect  on  the  results  which 
follow.  In  all  equilibrium  catalytic  wall  (ECW)  calculations, 
the  complete  multicomponent  equilibrium  wall  boundary  condi¬ 
tion  including  mass  injection  of  an  inert  specie  was  used.  It 
is  seen  that  helium  injection  had  a  very  small  influence  on 
the  heat  flux,  boundary- layer  thickness,  and  shock-layer  thick¬ 
ness  due  to  the  small  injection  rate.  Argon  had  a  larger  in¬ 
fluence  with  a  slight  thickening  of  both  boundary  and  shock 
layers  due  to  the  mass  addition. 

For  purposes  of  comparison  it  is  instructive  to  consider 
the  stagnation  point  heat  flux  for  the  thin  viscous  shock 
layer  relative  to  the  nonequilibrium  boundary- layer  using 
inviscid  nonequilibrium  outer  edge  conditions.  Such  has  been 
done  in  Table  VII  where  corresponding  wall  conditions  have 
been  matched,  e.g.,  boundary  layer  with  noncatalytic  wall  to 
thin  viscous  shock  layer  with  noncatalytic  wall  conditions. 

Also  shown  are  results  from  Goldberg  (Ref.  34),  who  employed 
a  viscous  shock  layer  model  in  complete  chemical  equilibrium, 
and  Cheng  (Ref .  28) ,  who  used  a  perfect  gas  viscous  shock 
layer  with  €  =  0.13,  M  and  a  highly  cooled  wall.  The 

increase  in  heat  transfer  at  the  stagnation  point  is  due 
mainly  to  vorticity  interaction  which  is  accounted  for  in  the 
thin  viscous  shock  layer  solution  but  neglected  in  the 
boundary- layer  solution.  This  vorticity  effect  manifests 
itself  through  a  finite  gradient  in  the  streamwise  velocity 
near  the  outer  edge  of  the  shock  layer.  Such  can  be  seen 
in  Fig.  23  which  shows  comparisons  of  the  thin  viscous  shock 
layer  velocity  and  temperature  profiles  with  respect  to 
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corresponding  nonequilibrium  boundary- layer  profiles  using 
nonequilibrium  inviscid  outer  edge  conditions.  The  fuller 
characteristic  of  the  boundary- layer  velocity  profile  is  a 
direct  consequence  of  imposing  a  zero  gradient  (no  external 
vorticity)  on  the  velocity  profile  at  the  outer  edge.  Al¬ 
though  Fig.  23  is  for  a  noncatalytic  wall  without  mass  injec¬ 
tion,  the  results  are  typical  for  all  cases,  even  including 
mass  injection.  Figure  24  shows  the  specie  distributions  at 
the  stagnation  point  from  the  viscous  shock  layer  analysis; 
note  that  both  wall  conditions  (with  and  without  mass  injec¬ 
tion)  are  presented.  The  important  fact  to  be  gained  from 
these  profiles  is  that  the  flow  field  remains  essentially 
frozen  in  all  cases.  There  are  no  regions  of  extreme  chemical 
production  or  depletion  due  to  nonequilibrium  phenomena,  and 
hence  multicomponent  diffusion  of  species  is  the  predominant 
mass  transport  mechanism.  Note,  however,  that  for  the  case 
of  an  equilibrium  catalytic  wall,  both  with  and  without  mass 
injection,  there  is  a  region  very  near  the  surface  which  is 
highly  influenced  by  the  wall  catalyticity  condition,  which 
is  to  be  expected. 

Much  of  the  foregoing  has  been  concerned  with  trying  to 
establish  whether  a  boundary-layer  analysis  is  valid  for  the 
present  body  and  set  of  conditions  and,  if  so,  what  the  cor¬ 
rect  outer  edge  conditions  should  be.  At  this  point  it  is 
well  to  summarize  our  findings  since  this  question  is  of 
fundamental  importance.  With  respect  to  the  outer  edge  con¬ 
ditions  for  the  boundary- layer  analysis,  it  was  shown  that 
the  stagnation  point  heat  flux  was  highly  dependent  on 
whether  the  outer  edge  was  taken  to  be  in  chemical  equilib¬ 
rium  or  nonequilibrium;  furthermore,  the  wall  catalyticity 
(equilibrium  catalytic  or  noncatalytic)  was  of  equal  impor¬ 
tance  since  it  controlled  the  energy  transfer  due  to  atom 
recombination  on  the  body  surface.  Based  on  this  information 
alone,  one  cannot  determine  the  more  correct  outer  edge  con¬ 
dition  and,  hence,  some  other  source  must  be  consulted.  The 
thin  viscous  shock  layer  provided  us  with  such  a  tool .  For 
the  body  and  set  of  conditions  under  examination,  it  was 
shown  that  a  thin  viscous  shock  layer  analysis  was  applicable. 
Comparison  of  stagnation  point  heat  fluxes  (shown  in  Table  II) 
revealed  that  the  nonequilibrium  outer  edge  boundary- layer 
conditions  were  in  good  agreement  with  the  thin  viscous  shock 
layer  for  both  wall  catalyticit ies  with  the  thin  viscous 
shock  layer  having  the  higher  heat  flux  value  due  to  the  in¬ 
clusion  of  vorticity  interaction  which  was  not  considered  in 
the  boundary-layer  treatment.  Furthermore,  it  was  shown  that 
the  boundary- layer  thickness  with  nonequilibrium  outer  edge 
conditions  was  about  65  percent  of  the  inviscid  shock  layer 
thickness,  i.e.,  <VAinv^sc^ci  ~  0.65,  as  well  as  being  thinner 
than  the  thin  viscous  shock  layer  thickness,  6/ApygL  ~  0.92. 
Based  on  these  results  it  is  felt  that  the  boundary-layer 
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analysis  which  used  the  nonequilibrium  outer  edge  conditions 
is  indeed  applicable  to  the  present  problem.  We  believe 
these  findings  and  conclusions  to  be  of  considerable  impor¬ 
tance  and  interest  since,  to  our  knowledge,  prior  investiga¬ 
tions  have  used  equilibrium  outer  edge  conditions  exclusively 
with  no  consideration  given  to  their  validity  on  physical 
grounds . 

The  above  discussion  was  concentrated  solely  on  the 
stagnation  point  where  it  was  shown  that  a  boundary-layer 
analysis  was  indeed  applicable  to  the  present  problem.  Hence, 
solutions  have  been  obtained  over  the  complete  body  using  the 
nonsimilar,  nonequilibrium  boundary- layer  program  described 
earlier.  These  results  will  now  be  presented. 

4.2.2  Downstream  Region 


One  of  our  primary  objectives  in  this  work  concerned  the 
effects  of  inert  gas  injection  on  various  boundary- layer 
parameters.  Hence,  a  measure  of  the  inert  gas  concentration 
at  the  wall  is  of  interest  since  it  essentially  controls  the 
inert  gas  influence  on  the  boundary  layer.  Figure  25  shows 
the  wall  concentration  of  the  air  species  defined  by  Zw  =  1  - 
cw  inert  gas5  and  thus  more  inert  gas  present  at  the  wall, 
the  smaller  the  value  for  Zw,  The  important  point  to  note 
here  is  the  large  influence  of  wall  catalyticity  on  inert  gas 
concentration  over  the  downstream  region  of  the  body.  Near 
the  end  of  the  body  the  noncatalytic  wall  condition  gives  a 
40-  to  85-percent  increase  in  the  inert  gas  concentration  at 
the  wall  as  compared  to  the  equilibrium  catalytic  wall  condi¬ 
tion.  This  is  both  surprising  and  enlightening. 

One  would  expect  the  mass  injection  process  to  have  a 
thickening  effect  on  the  boundary- layer  displacement  thick¬ 
ness.  Such  is  true  as  shown  in  Fig.  26.  Note  first  that  for 
the  case  of  no  mass  injection,  the  effect  of  wall  catalyticity 
is  negligible.  Such  is  not  true  with  inert  gas  injection; 
here  the  noncatalytic  wall  condition  has  a  much  larger  dis¬ 
placement  effect  than  does  the  equilibrium  catalytic  wall  in 
the  downstream  region.  Considering  the  way  that  viscous  inter¬ 
action  involves  the  displacement  thickness,  one  should  expect 
larger  viscous  interaction  effects  for  the  noncatalytic  wall 
condition  than  for  the  equilibrium  catalytic  wall  with  mass 
injection. 

Turning  now  to  the  heat-transfer  results,  Figs.  27  and 
28  show  the  heat  flux  distributions  for  all  conditions  over 
the  entire  body.  It  is  seen  that  the  noncatalytic  wall  re¬ 
duces  the  heat  flux  as  expected  and  by  a  substantial  amount 
in  the  case  of  argon  injection.  It  is  interesting  to  note 
that  the  difference  between  the  equilibrium  catalytic  wall  and 
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noncatalytic  wall  results  remains  essentially  constant  over 
the  entire  body  for  the  no-mass-injection  condition.  For 
cases  with  mass  injection,  the  type  of  plot  shown  in  Fig.  29 
is  instructive  since  it  shows  directly  the  reduction  in  heat 
flux  due  to  mass  addition  relative  to  the  no-injection  result 
for  a  given  wall  catalyticity .  Here  the  effectiveness  of 
argon  injection  with  a  noncatalytic  wall  becomes  immediately 
apparent.  Hence,  it  is  clear  that  it  is  highly  desirable  to 
have  noncatalytic  surfaces  if  one  desires  to  substantially 
reduce  the  wall  heat  transfer  by  normal  injection  of  an  inert 
gas  (transpiration  cooling)  into  a  nonequilibrium  boundary 
layer . 

It  is  also  of  interest  to  examine  the  skin-friction 
distribution  shown  in  Fig.  30.  For  the  case  of  no  mass  in¬ 
jection,  there  is  no  influence  due  to  wall  catalyticity  on 
the  skin  friction.  However,  when  mass  injection  is  consid¬ 
ered,  the  noncatalytic  wall  reduces  the  skin  friction  as  com¬ 
pared  with  the  equilibrium  catalytic  wall.  Since  the  skin- 
friction  drag  results  from  integrating  the  wall  shear  stress 
distribution  over  the  body,  the  possibility  of  reducing  such 
drag  by  injection  of  an  inert  gas  is  well-known.  Table  VIII 
shows  this  reduction  in  skin-friction  drag  with  inert  gas 
injection;  as  expected,  the  noncatalytic  wall  condition  re¬ 
sults  in  the  lower  drag  for  both  gases. 

In  summary,  the  above  discussion  has  pointed  out  the 
possibility  of  reducing  both  heat  transfer  and  skin  friction 
by  injection  of  an  inert  gas  (argon  or  helium  in  this  case) . 

In  addition,  a  further  reduction  in  both  heat  transfer  and 
skin  friction  is  possible  with  such  inert  gas  injection  by 
using  a  noncatalytic  surface  as  opposed  to  an  equilibrium 
catalytic  surface.  However,  such  a  choice  will  increase  the 
resulting  boundary- layer  displacement  thickness  which  could 
result  in  larger  viscous  interaction  effects.  Hence,  some 
sort  of  a  "tradeoff”  may  prove  desirable  in  future  applications. 


4.3  COMPARISON  OF  BINARY  GAS  AND  NONEQUILIBRIUM  FLOW  FIELD  RESULTS 


A  comparison  of  the  binary  gas  and  nonequilibrium  ioniz¬ 
ing  air  boundary- layer  results  is  shown  in  Fig.  31.  For  pur¬ 
poses  of  this  comparison,  binary  gas  calculations  were  made 
without  transverse  curvature  and  displacement  effects. 

The  comparison  shows  good  agreement  for  all  the  quanti¬ 
ties  plotted  in  Fig.  31  for  the  case  of  no  mass  injection. 
Smaller  differences  in  the  calculated  results  were  found  for 
the  case  of  helium  injection  than  for  argon  injection.  The 
effects  of  the  higher  injection  rate  of  argon  on  the 
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displacement  thickness  comparison  can  be  seen  in  Fig.  31a. 

Some  of  the  difference  in  the  calculated  heat- transfer  rates 
is  removed  when  the  slight  variation  in  free-stream  condi¬ 
tions  is  included.  This  can  be  seen  by  comparison  of  the 
differences  in  calculated  dimensional  heat-transfer  rate  in 
Fig.  31c  with  the  dimensionless  heat- transfer  rate  4/p^U^ 
shown  in  Fig.  31e . 

The  relatively  good  agreement  between  these  binary  gas 
and  ionizing  air  boundary-layer  results  should  support  the 
confidence  in  the  binary  gas  results  including  transverse 
curvature  and  displacement.  The  scope  of  these  calculations 
did  not  include  ionizing  air  boundary-layer  calculations 
with  higher-order  boundary-layer  effects.  The  reacting  gas 
boundary-layer  calculations  were  made  to  assess  the  effects  of 
nonequilibrium  chemistry  without  higher-order  boundary- layer 
effects  and  thus  would  be  useful  in  interpreting  the  binary 
gas  boundary-layer  calculations.  The  comparisons  in  Fig.  31 
show  a  measurable  nonequilibrium  effect  on  the  heat  transfer 
with  argon  injection  near  the  end  of  the  conical  afterbody. 
Some  of  the  difference  between  binary  gas  and  ionizing  air 
calculated  heat- transfer  results  might  be  reduced  if  trans¬ 
verse  curvature  and  displacement  were  included  in  the  ioniz¬ 
ing  air  model.  Therefore,  the  effects  of  nonequilibrium 
chemistry  cannot  be  predicted  accurately  from  these  calcula¬ 
tions  and  comparisons. 


4.4  EFFECTS  OF  WALL  TEMPERATURE 

The  effects  of  Avail  temperature  on  the  stagnation  point 
heat-transfer  rate  for  the  binary  gas  and  ionizing  air  boundary- 
layer  models  without  mass  transfer  are  shown  in  Fig.  32.  Since 
the  experimental  no-injection  wall  temperature  may  not  have 
been  the  assumed  value  of  540°R  (the  same  as  the  wall  tempera¬ 
ture  conditions  with  injection),  it  was  considered  desirable 
to  investigate  the  effects  of  wall  temperature  variation.  The 
change  in  heat- transfer  rate  with  wall  temperature  was  small 
for  540  <  Tw  <  700°R  and  the  rate  of  change  was  approximately 
constant  for  "500  <  Tw  <  1500°R  where  the  upper  limit  was  the 
maximum  value  considered.  With  the  data  shown  in  Fig.  32,  it 
is  possible  to  correct  the  ratio  4inj/4no  inj  shown  in  Figs. 

22  and  29  for  the  effect  of  wall  temperature.  It  is  interest¬ 
ing  to  note  that  for  the  range  of  wall  temperatures  considered 
the  difference  between  equilibrium  catalytic  wall  (ECW)  and 
noncatalytic  wall  (NCW)  heat  transfer  remained  essentially 
constant,  and  similarly  for  Tw  >  800°R  the  difference  between 
the  nonreacting  binary  gas  and  ionizing  air  with  NCW  was  in¬ 
sensitive  to  wall  temperature. 
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SECTION  V 

CONCLUDING  REMARKS 


Binary  gas  laminar  boundary- layer  calculations  were  made 
including  mass  transfer,  transverse  curvature,  and  displace¬ 
ment  effects.  Nonequilibrium  ionizing  air  boundary  layer  and 
stagnation  point  thin  viscous  shock  layer  calculations  were 
made  including  mass  transfer  and  catalytic  wall  effects.  The 
comparison  of  binary  gas  and  ionizing  air  boundary- layer 
results  including  mass  transfer  but  excluding  higher-order 
boundary- layer  effects  showed  reasonably  good  agreement  except 
in  the  case  of  argon  injection  and  supported  the  use  of  the 
binary  gas  results  including  transverse  curvature  and  dis¬ 
placement  effects.  The  data  were  presented  to  aid  in  making 
comparisons  between  the  calculated  results  and  experimental 
data.  Final  assessment  of  the  applicability  of  the  theoret¬ 
ical  models  used  and  the  numerical  results  obtained  must  await 
such  comparison  with  experimental  data  currently  unavailable. 
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Fig.  19  Dimensionless  Heat-Transfer  Rate  q/p^U^5  to  the  Cone  with  Helium  Injection 
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Fig.  20  Dimensionless  Heat-Transfer  rate  q/q(0)  to  the  Cone  with  Helium  Injection 
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Fig.  23  Velocity  and  Temperature  Profiles  at  the  Stagnation 
Point  under  Nonequilibrium  Conditions 
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y 

b.  Argon  Injection 
Fig.  24  Continued 
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Fig.  31 


a.  Displacement  Thickness 

Comparison  of  Binary  Gas  and  Ionizing  Air  Boundary-Layer  Results  without  Transverse 
Curvature  or  Displacement  Effects:  =  13.2 ,6.  =  7.5  deg,  T„  =  540°R 
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c.  Dimensional  Heat-Transfer 


Fig.  31  Continued 
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e.  Dimensionless  Heat-Transf 


Fig.  31  Concluded 
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Fig.  32  effects  of  Wall  Temperature  on  No-Injection  Stagnation  Heat-Transfer  Rate 
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TABLE  I 

FREE-STREAM  AND  STAGNATION  CONDITIONS* 


* 

Binary  Gas 
Calculations 

Nonequilibrium 

Calculations 

M 

00 

13.2 

— 

Poo’ 

lbf/ft2 

0.1944 

0.1944 

u  , 

00 

ft/sec 

10,680 

10,850 

T  , 

CO 

°R 

290.2 

274 

Poo’ 

slugs/f t^ 

3.8998  x  10"7 

4.0716  x  10 

p  U  ,  lbf-sec/ft^ 
r  00  00  7 

0.004165 

0.004418 

Po’ 

lbf/ft2 

41.647 

45.886 

To’ 

°R 

8204 

8264.2 

Tw’ 

°R 

540 

540 

Ho’ 

ft2/sec^ 

5.9047  x  107 

6.6011  x  10 

V 

1.313 

_  _ 

*Free-stream  chemical  composition  data,  which  are 
equivalent  to  the  frozen  normal  shock  specie  concentra¬ 
tions,  are  given  in  Tables  II,  III,  and  IV. 
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TABLE  II 

NONEQUILIBRIUM  STAGNATION  POINT  CONDITIONS  WITH  NO  MASS  INJECHON 


Boundary-Layer  Solution 

NEQ  In vise id  Outer  Edge 

Boundary-Layer  Solution 

EQ  Invlscld  Outer  Edge 

Thin  Viscous  Shock  Layer 

Frozen  Normal  Shock  Crossing 

Wall 

(NEQ, FCW) 

Wall 

(NEQ ,  NCW  ) 

Outer 

Edge (NEQ) 

Wall 

(NEQ, FCW) 

Wall 

(NXQ.NCW) 

Outer 

Edge(EQ) 

Wall 

(NEQ, FCW) 

wall 

(NEQ, NCW) 

Behind  Shock 
(Frozen) 

Q 

0.0 

2.3552-02 

2.3186-02 

0.0 

2.3339-01 

2 . 3146-01 

0.0 

2.0478-02 

2.0432-02 

°2 

2 . 3456-01 

1.7634-01 

1.7701-01 

2.3456-01 

3.2498-04 

2.1200-04 

2.3456-01 

1.7996-01 

1.8000-01 

MO 

0.0 

6.5063-02 

6.4895-02 

0.0 

8.7191-04 

5.3410-03 

0.0 

6.4441-02 

6.4455-02 

H 

0.0 

1.3513-05 

8.2880-05 

0.0 

3.3963-02 

3.5750-02 

0.0 

7 . 9708-06 

1.9295-09 

NO* 

0.0 

5.1212-09 

3.5248-10 

0.0 

2.6158-05 

6.4000-05 

0.0 

3.6449-08 

3.6459-08 

N2 

7.6544-01 

7.3504-01 

7.3482-01 

7.6544-01 

7 . 3143-01 

7.2717-01 

7.6544-01 

7.3512-01 

7.3537-01 

[A]o 

(Btu/f t  -sec) 

55.94 

49.61 

101.63 

51.22 

66.68 

59.63 

6  (at 

£'  ■=  0.995) 

0.0695 

0.0696 

0.0812 

0.0819 

A 

0.1091 

0. 1095 

0.0753 

0.0758 

p^Clbf/ft2) 

45.886 

41.647 

41.471* 

(°R) 

8264.2 

8204.0 

8069.2** 

NEQ  (Nonequillbrium)  *p_  (lbf/ft2) 

KQ  (Equilibrium) 

FCW  (Fully  Catalytic  Wall)  **T  (”») 

HCW  (Noncatalytic  Wall) 


TABLE  HI 


NONEQUILIBRIUM  STAGNATION  POINT  CONDITIONS  WTH  ARGON  INJECTION  AT 
.Pwvw  =  1-182  x  10-4  lbf-soc/ft^ 


Boundary-Layer  Solution 

NEQ  Inviscid  Outer  Edge 

Thin  Viscous  Shock  Layer 

Frozen  Normal  Shock  Crossing 

Wall 

(NEQ, ECW) 

Wall 

(NEQ.NCW) 

Outer 

Edge (NEQ) 

Wall 

(NEQ, ECW) 

Wall 

(NEQ,NCW) 

Behind  Shock 
(Frozen) 

0 

°2 

2.0438-40 

1.9646-01 

2.0341-02 

1.4039-01 

2 . 3186-02 
1.7701-01 

2.0774-40 

1.8219-01 

1.8196-02 

1.4905-01 

2.0432-02 

1.8000-01 

NO 

N 

NO+ 

Argon  (Inj) 

N2 

7.4851-16 

2.0710-79 

8.1040-86 

1.6157-01 

6.4197-01 

5.4120-02 

1.2308-05 

4.0881-09 

1.7168-01 

6.1346-01 

6.4895-02 

8.2880-05 

3.5248-10 

0.0 

7.3482-01 

7.3557-16 

2.2299-79 

8.4785-86 

1.4936-01 

6.6845-01 

5.5391-02 

7.5081-06 

3.2054-08 

1.4442-01 

6.3293-01 

6.4455-02 

1.9295-09 

3.6459-08 

0.0 

7.3537-01 

["iJo 

(Btu/f t2-sec) 

51.29 

45.56 

61.82 

55.20 

6  (at 

f'  =  0.995) 

0.0720 

0.0720 

A 

0.1091 

0.0770 

0.0773 

p^( lbf/ft2) 

45.886 

41.471* 

(°R) 

8264.2 

8069.2** 

ECW  (Equilibrium  Catalytic  Wall) 


*pe  (lbf/ft2) 
s 

**T  (°R) 

s 
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TABLE  IV 

NONEQUILIBRIUM  STAGNATION  POINT  CONDITIONS  WITH  HELIUM  INJECTION  AT 
pwv w  ,  3.080  X  10~5  Ibf-xee/ft3 


Boundary-Layer  Solution 

NEQ  Inviscid  Outer  Edge 

Thin  Viscous  Shock  Layer 

Frozen  Normal  Shock  Crossing 

Wall 

(NEQ,ECW) 

Wall 

(NEQ,NCW) 

Outer 

Edge (NEQ) 

Wall 

(NEQ,ECW) 

Wall 

(NEQ,NCW) 

Behind  Shock 
(Frozen) 

0 

°2 

2.3607-40 

2.2216-01 

2.3292-02 

1.7118-01 

2.3186-02 

1.7701-01 

2.3645-40 

2.1250-01 

2.0492-02 

1.7690-01 

2.0432-02 

1.8000-01 

NO 

N 

NO+ 

Helium  (Inj) 

N2 

8.6475-16 

2.4428-79 

9.4592-86 

2.0138-02 

7.5770-01 

6.3111-02 

1.3038-05 

4.4639-09 

2.0638-02 

7.2177-01 

6.4895-02 

8.2880-05 

3.5248-10 

0.0 

7.3482-01 

8.5652-16 

2.5337-79 

9.6420-86 

1.0434-02 

7.7706-01 

6.3385-02 

7.0991-06 

3.5020-08 

1.5673-02 

7.2355-01 

6.4455-02 

1.9295-09 

3.6459-08 

0.0 

7.3537-01 

(Btu/it2-see) 

55.81 

49.46 

65.41 

58.76 

6  (at 

f  '  -  0.995) 

0.0651 

0.0652 

A 

0.1091 

0.075? 

0.0759 

p^(lbf/ft2) 

45.886 

41.471* 

(°R) 

8264.2 

8069.2** 

*pc  (lbf/ft2) 

S 


**T  (°R) 
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TABLE  V 

CHEMICAL  REACTIONS  AND  RATES 


Reaction 

Catalytic  Body 
<M) 

Cf 

af 

e  f/k 
0,1 

(°K) 

Cb 

ab 

e  ./k 
o ,  b 

<°K) 

02  +  °2  -  20  +  02 

2.3  x  1019 

-1 

59,400 

1.9  x  1016 

-1/2 

0 

02  +  0  2  20  +  0 

8.5  x  1019 

-1 

59,400 

7.1  x  1016 

-1/2 

0 

02  +  M  2  20  +  M 

N,N2,NO, Inert 

3.0  x  1018 

-1 

59,400 

2.5  x  1015 

-1/2 

0 

N2  +  N2  2  2N  +  N2 

3.8  x  1019 

-1 

113,200 

2.0  x  1018 

-1 

0 

N2  +  N  2  2N  +  N 

1.3  x  102° 

-1 

113,200 

7.0  x  1018 

-1 

0 

N2  +  M  2  2N  +  M 

0,02 ,NO, Inert 

1.9  x  1019 

-1 

113,200 

1.0  x  1018 

-1 

0 

NO  +  M  2  N+O+M 

0,02>N,N2, Inert 

2.4  x  1017 

-1/2 

75 , 500 

6.0  x  1016 

-1/2 

0 

NO  +  O  2  02  +  N 

4.3  x  107 

3/2 

19,100 

1.8  x  108 

3/2 

3020 

N2  +  O  2  NO  +  N 

6.8  x  1013 

0 

37,750 

1.5  x  10Ad 

0 

0 

n2  +  °2  *  2N0 

14 

2.0  x  10  * 

0 

61,600 

1.0  x  1013 

0 

40,000 

N  +  O  2  NO+  +  e” 

1.3  x  108 

1 

31,900 

2.0  x  1019 

-1 

0 

Note:  Reaction  rates  are  from  Bortner,  -M.  H.  "Chemical  Kinetics  in  a  Re-Entry 

Flow  Field."  GE  TIS  R63SD63,  August  1963. 


3 

kf  =  Cf  T  exp  (-eQ  f /kT) ,  cm  /mole-sec  kb  "  Cb  T  exp  (-eQ  b/kT) 

T  in  °K 
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TABLE  VI 

COMPARISON  OF  FULLY  CATALYTIC  WALL  AND  EQUILIBRIUM  CATALYTIC 
WALL  SPECIE  MASS  FRACTIONS  UNDER  NONEQUILIBRIUM  CONDITIONS 


Specie 

Fully  Catalytic 
Wall  (FCW) 

Equilibrium  Catalytic 
Wall  (ECW) 

0 

0.0 

2.4038-40 

°2 

2.3456-01 

2.3451-01 

NO 

0.0 

8.9306-16 

N 

0.0 

2.4345-79 

NO+ 

0.0 

9.5297-86 

N2 

7.6544-01 

7.6549-01 

Note:  Above  conditions  based  on  thin  viscous  shock 

layer  analysis  at  the  stagnation  point  without 
mass  transfer. 
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TABLE  VII 

STAGNATION  HEAT-TRANSFER  RATE  UNDER  NONEQUILIBRIUM  CONDITIONS 


c^,  Thin  Viscous-Shock  Layer 


Boundary-Layer 


No  Injection 

NEQ,  FCW  1.192 

NEQ ,  NCW  1 . 202 

Goldberg  (Ref.  34)  «  1.19 

Cheng  (Ref.  28)  *1,13 

Argon  Injection 

NEQ,  ECW  1.205 

NEQ,  NCW  1.212 

Helium  Injection 

NEQ,  ECW  1.172 

NEQ,  NCW  1.188 
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TABLE  VIII 

MASS  TRANSFER  EFFECTS  ON  DRAG  REDUCTION 


Gas  Model 

Displacement 

Iteration 

TVC 

In jectant 

Noncatalytic 

Wall 

Equilibrium 

Wall 

Reacting 

0 

No 

Argon 

0.488 

0.582 

Reacting 

0 

No 

Helium 

0.671 

0.759 

Nonreacting 

0 

No 

Argon 

0.4494 

— 

Nonreacting 

0 

No 

Helium 

0.6501 

— 

Nonreacting 

0 

Yes 

Argon 

0.4502 

— - 

Nonreacting 

0 

Yes 

Helium 

0.6578 

— 

Nonreacting 

5 

Yes 

Argon 

0.5299 

— 

Nonreacting 

3 

Yes 

Helium 

0.7395 

— 

♦Values  Given  for  Ratio  CL.  /CL 

i  '  Uf 
inj  no  inj 
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APPENDIX  III 

FLUID  PROPERTIES  USED  IN  THE  BINARY  GAS 
BOUNDARY-LAYER  PROGRAM 


To  permit  computations  for  which  temperatures  in  the 
boundary  layer  exceed  6300°R,  it  was  necessary  to  extend  the 
range  of  validity  of  the  fluid  properties  in  the  Jaffe,  Lind, 
and  Smith  boundary-layer  computer  program  (Ref.  5).  The 
basic  properties  required  are  cp  and  p.  for  each  individual 
species  and  Djj  for  the  binary  pairs  of  interest.  Other 
needed  speciesand  mixture  properties  are  calculated  from 
these  basic  properties  in  the  manner  described  in  Ref.  5. 

The  fluid  properties  have  been  extended  to  12,600°R  for  the 
free-stream  species  air  and  the  injectant  species  helium  and 
argon. 

The  constant  pressure  specific  heat  for  air  from  0  to 
12,600°R  was  taken  from  Dommett’s  data  for  undissociated  air 
(Ref.  38),  and  was  fit  with  two  fifth-degree  polynominals, 
one  from  0  to  2000°R  and  another  from  2000  to  12,600°R.  The 
form  of  the  polynominals  and  the  values  of  the  coefficients 
are  given  in  Table  III-l. 

For  helium  and  argon,  only  the  contribution  of  the  trans¬ 
lational  mode  was  considered  in  determining  the  specific  heats. 
This  results  in  a  constant  value  of  cp  for  each  of  the  two 
species  of 

cp  =  (5/2)  R 

This  does  not  involve  any  change  from  the  data  in  Ref.  5,  only 
an  extension  of  the  temperature  range  to  which  it  is  applied. 

The  viscosity  of  air  from  90  to  720°R  was  computed  from 
Sutherland’s  equation,  and  from  900  to  12,600°R  the  data  from 
Peng  and  Pindroh  (Ref.  39)  at  a  density  of  0.1  amagat  were 
used.  These  two  sets  of  data  were  smoothed  together  and  fit 
with  a  fifth-degree  polynominal.  The  form  of  the  polynominal 
and  ,the  coefficients  are  given  in  Table  III-2. 

The  viscosities  of  helium  and  argon  as  given  in  Ref.  5 
were  used  up  to  3000°R;  from  3000  to  12,600°R  the  data  of 
Amdur  and  Mason  (Ref.  40)  were  used.  These  data  were  also 
fit  with  fifth-degree  polynominals,  and  the  coefficients  are 
given  in  Table  III-2. 

The  binary  diffusion  coefficients  for  helium  and  argon 
in  air  were  computed  using  the  Chapman-Enskog  relation  and  the 
Lenard-Jones  (6-12)  intermolecular  potential  function.  The 
force  constants  for  the  individual  species  were  obtained  from 
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Hirschf elder ,  Curtiss,  and  Bird  (Ref.  41).  The  resulting 
diffusion  coefficients  were  fitted  with  fifth-degree  poly- 
nominals  as  given  in  Table  III-3. 

Plots  of  the  various  fluid  properties  just  described  are 
given  in  Figs.  III-l,  III-2,  and  III-3.  These  plots  are  from 
the  values  of  the  properties  as  computed  by  the  appropriate 
fifth-degree  polynominals  and  not  the  original  data. 
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TABLE  lll-l 

CONSTANT  PRESSURE  SPECIFIC  HEAT  OF  AIR 
cp  “  A  +  BT  +  CT2  +  DT3  +  ET4  +  FT5  (f t2/sec2-°R) 


Temperature 

0  -  2000°R 

2000  -  12, 600°R 

A 

6.0371797  x  103 

5,9028  x  103 

B 

-9.4509125  x  10-4 

3.77072  x  10_1 

C 

-7.3022675  x  10-4 

9,64649  x  10-5 

D 

1.73022675  x  10"6 

-3.53769  x  10“8 

E 

-9.7657438  x  10"10 

3.48567  x  10"12 

F 

1.7465179  x  10-13 

-1.11502  x  10"16 

TABLE  111*2 

VISCOSITY  OF  AIR,  HELIUM,  AND  ARGON 
|X  x  107  -  A  +  BT  +  CT2  +  DT3  +  ET4  +  FT5  (lbf-sec/f t2) 


Gas 

Air 

Helium 

Argon 

A 

1.48066  x  10  1 

7.2044  x  10-1 

2.63154  x  10-1 

B 

6.95936  x  10“3 

7.06794  x  10“3 

8.61381  x  10-3 

C 

-1.49079  x  10-6 

-1.5363  x  10'6 

-1.84422  x  10-6 

D 

2.3759  x  10"10 

2.80513  x  10-1° 

3.16427  x  10_1° 

E 

-1.78242  x  10-14 

-2.28363  x  10'14 

-2.47897  x  10-14 

F 

5.0725  x  10"19 

6.74097  x  10-19 

7.10697  x  10-19 

TABLE  111*3 

BINARY  DIFFUSION  COEFFICIENTS 

pDi;j  »  A  +  BT  +  CT2  +  DT3  +  ET4  +  FT5  (lbf/sec) 


;  Species 

Helium-Air 

Argon-Air 

A 

-1.98803  x  10"1 

-6.39025  x  10-2 

B 

2.31693  x  10“3 

6.67803  x  10-4 

C 

2.60637  x  10"6 

1.26081  x  10-6 

D 

-4.74411  x  10"11 

-1.02832  x  10"10 

E 

-1.00312  x  10“14 

7.39182  x  10"15 

F 

6.79428  x  10“19 

-2.18881  x  10"19 
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Fig.  111-2  Viscosity  of  Air,  Helium,  and  Argon 
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APPENDIX  IV 

i 

SURFACE  HEAT-TRANSFER  RELATIONS 


For  the  case  of  nonreacting  gas  flow  over  a  nonadiabatic 
surface  with  no  mass  injection,  the  resulting  wall  heat  flux 
is  given  by  the  well-known  Fourier  law  of  heat  conduction  as 


(IV-1) 


However,  when  the  flow  is  chemically  reacting  or  mass  transfer 
is  allowed,  the  classical  definition  above  must  be  modified  to 
include  these  effects.  This  appendix  will-  be  devoted  to  de¬ 
riving  a  general  surface  heat- transfer  expression  which  allows 
all  of  these  effects  to  occur  simultaneously  and  hence  is  valid 
for  both  reacting  and  nonreacting  flows  either  with  or  without 
mass  injection. 


Following  Chung  (Ref.  21),  the  surface  is  assumed  to  be  a 
rigid  semipermeable  membrane  which  does  not  allow  the  air  species 
(denoted  by  subscript  f)  to  flow  into  the  porous  surface  but  does 
allow  the  injected  species  (denoted  by  subscript  i)  to  flow  into 
the  boundary-layer  or  shock  layer.  At  the  surface,  energy  is 
transported  from  the  boundary-layer  or  shock  layer  to  the  surface 

bT 


by  conduction  and  diffusion  at  the  rate  ^k  ^  -  2  jaha)w 


as  shown 


in  the  sketch  below. 


37  - 1 

At  the  same  time  energy  is  transported  away 


from  the  surface-gas  interface  by  the  normal  mass  flow  (pv)w  at 
the  rate  (pv)w  g  Caw  haw  and  toward  the  solid-gas  interface  at 

the  rate  (pv)w  2  (Cpw  hpw)_,  The  subscript  (-)  denotes  the  in¬ 


terior  of  the  solid,  whereas  the  subscript  (P)  denotes  the  injec¬ 
ted  species  i.  The  subscript  (a)  denotes  the  air  species  f  as 
well  as  the  injected  species  i.  In  order  to  clarify  this  nota¬ 
tion  let  us  consider  a  six  species  air  model  (0,  O2,  N,  N2,  NO, 
N0+)  with  injection  of  helium  (He).  For  this  example 


f  =  0,  02,  N,  N2,  NO,  N0+ 
i  =  He 

a  -  0,  02,  N,  N2,  NO,  N0+, 


He  (a 


P  =  He  (p  =  i) 


f  +  i) 


where  (Cpw) _  =1.0  since  this  quantity  is  evaluated  in  the  in¬ 
terior  of  the  solid  and  only  one  species  is  being  injected.  In 


79 


AEDC-TR-68-214 


addition,  2  (Cpw)_  =1.0  and  2  Caw  =  1.0  are  general  restric¬ 
tions  whicfi  must  be  imposed  in  all  cases. 


Gas-Solid 
Interface 
Treated  as 
Semipermeable 
Membrane 


Based  on  the  above  sketch,  the  net  rate  of  heat  transfer 
into  the  solid  is 


(IV-2) 


where  a  (+)  value  for  qw  denotes  heat  transfer  into  the  interior 
of  the  solid.  Note  that  for  the  case  of  a  nonreacting  single 
species  gas  flow  with  no  mass  injection  Eq.  (IV-2)  reduces  to 
Eq.  (IV-1) . 


In  order  to  reduce  Eq.  (IV-2)  to  a  simpler  form,  it  is 
necessary  to  examine  a  mass  transfer  balance  at  the  gas-solid 
interface  as  shown  in  the  sketch  below. 


Gas-Solid 
Interface  Treated 
as  Semipermeable 
Membrane 
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The  component  a  is  transported  from  the  gas  to  the  solid  by 
diffusion  at  the  rate  -jaw.  At  the  same  time  the  component 
a  is  transported  away  from  the  interface  by  the  normal  mass 
flow  (pv)w  at  the  rate  (pv)w  Caw  while  the  component  a  may 
possibly  be  transported  toward  the  interface  from  within  the 
solid  at  the  rate  (pv)w  (Ca*  )_  where  if  P*  ^  a,  then 
(Cfi*w) -  =  0 •  If  P*  =  a,  then  (Cp*w)_  =  (Caw)_  where  the  (-) 
subscript  is  to  be  interpreted  as  discussed  previously.  Sum¬ 
ming  these  components  gives  the  net  rate  of  production . of  the 
component  a  by  surface  chemical  reactions  (denoted  by  Sa)  as 


Sa  =  Jaw  +  <pv)w 


"  <cp*w)_] 


(IV-3) 


For  the  case  of  a  noncatalytic .surface  which  allows  no  atom 
recombination  on  the  surface,  Sa  =  0  for  each  component  a. 

For  cases  in  which  heterogeneous  chemical  reactions  are  allowed 
to  occur  at  the  gas-solid  interface,  i.e.,  catalytic  surface 
recombination  of  atoms,  a  detailed  treatment  of  the  surface 
chemical  kinetics  is  required  in  order  to  evaluate  Sa.  Since 
the  detailed  theoretical  treatment  of  the  surface  chemical 
kinetics  is  even  more  complex  than  that  of  the  homogeneous 
ones,  it  will  not  be  considered  in  the  present  work.  The 
reader  is  referred  to  Chung  (Ref.  21)  for  more  details  on  this 
subject. 

Returning  to  Eq.  (IV-3),  note  that  conservation  of  mass 
requires  2  Sa  =  0  and  hence  2  =  0  since  2  CaTO  =1.0  and 

a  CL  w  a 

2  (Cp*w)_  =  1.0.  For  the  case  of  a  noncatalytic  wall  (Sa  =  0) , 


Eq.  (IV-3)  reduces  to 


Ja»  "  (Pv)w  [(cpV-  ‘  Caw] 


(IV-4) 


which  yields  the  well-known  noncatalytic  wall  boundary  condi¬ 
tion  jctw  =  0  for  cases  with  no  mass  injec'tion.  Since  jaw  « 

/dc  a\ 

l  — — )  Eq.  (IV-4)  relates  the  concentration  gradients  to  the 
\  oy  /w 

actual  concentrations  for  cases  with  mass  injection  through 
a  noncatalytic  surface.  Furthermore,  note  that 

^a^aw  “  jaw  ^w  +  ^pv^w  ^w  [^aw  - 

which  results  from  multiplying  Eq.  (IV-3)  by  the  static  enthalpy 
of  species  a  evaluated  at  the  wall.  Summing  the  above  relation 
over  all  species  yields 
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a 


Ja 


w  ~  ^-4  J 

a 


aw 


w 


+  <pv)w  IXw  [Caw  -  (Cp*»)J  (IV-5) 

and  recalling  that,  due  to  the  relationship  between  a ,  p* ,  and 
P  discussed  previously, 

E  <Vw>-  -  Ehp„  (Cp.)_  (IV-6) 

Q*  p 

Eq.  (IV-2)  becomes,  upon  substitution  of  Eqs.  (IV-5)  and  (IV-6), 


(IV-7) 

which  is  an  alternate  form  of  the  relation  for  the  net  rate 
of  heat  transfer  into,  the  interior  of  the  solid.  The  physical 
interpretation  of  Eqi  (IV-7)  is  clear:  the  right-hand  side  is 
the  difference  between  the  rate  of  heat  conduction  to  the  wall 
from  the  gas  flow  and  the  rate  of  heat  absorption  at  the  wall 
due  to  surface-catalyzed  reactions,  and  the  left-hand  side  is 
the  total  heat  transfer  into  the  interior  of  the  solid.  It  is 
interesting  to  note  that  the  heat  carried  to  the  surface  by 
diffusion  is  precisely  balanced  by  convection  of  heat  due  to 
mass  transfer  away  from  the  surface  and,  hence,  makes  no  net 
contribution  to  the  "sensible"  heat  balance. 

.  In  order  to  evaluate  Eq.  (IV-7),  the  net  rate  of  production 
Sa  must  be  known  for  each  component  a  which,  as  discussed  pre¬ 
viously,  is  not  easily  evaluated  for  cases  with  surface  cata- 
lyticity  effects.  However,  for  case  of  a  noncatalytic  surface, 
recall  that  Sa  =  0  for  each  component  a,  and  hence  Eq.  (IV-7) 
reduces  to  the  classical  form  of  the  Fourier  law  of  heat  con¬ 
duction,  Eq.  (IV-1) ,  irrespective .of  mass  injection  effects. 

This  is  an  important  point  since  Sa  =  0  for  each  a  in  the  case 
of  a  nonreacting  flow;  this  shows  that  one  should  properly 
compare  nonreacting  flow  heat  transfer  with  chemically  re¬ 
acting  noncatalytic  wall  heat  transfer  in  order  for  the  com¬ 
parison  to  have  physical  meaning.  Furthermore,  it  is  extremely 
important  to  note  that  mass  transfer  has  no  direct  effect  on 
the  expression  for  the  wall. heat  transfer,  Eq.  (IV-1),  which 
is  valid  for  the  condition  Sa  =  0  irrespective  of  mass  transfer. 

In  order  to  analyze  flows  with  wall  catalyticity  effects, 
either  Eq.  (IV-2)  or  Eq.  (IV-7)  must  be  used.  One  of  the  most 
common  assumptions  with  regard  to  a  catalytic  surface  is  that 
it  be  in  local  thermodynamic  equilibrium  at  the  given  wall 
temperature  and  pressure.  As  discussed  by  Moore  (Ref.  23)  such 
an  equilibrium  catalytic  wall  condition  requires  a  treatment  of 
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local  chemical  equilibrium  at  the  gas-solid  interface  including 
the  effects  of  diffusion  and  mass  injection.  Since  Moore  de¬ 
scribes  this  procedure  in  great  detail,  the  reader  is  referred 
to  his  work  for  a  complete  description. 

In  summary,  Eqs.  (IV-2)  and  (IV-3)  form  a  general  set  of 
heat  and  mass  transfer  relations  at  the  gas-solid  interface 
including  the  effects  of  mass  injection  through  a  porous  sur¬ 
face  considered  as  a  rigid  semipermeable  membrane  which  does 
not  allow  the  air  species  to  flow  into  the  porous  surface  but 
does  allow  the  injected  species  to  flow  into  the  boundary-layer 
or  shock  layer.  For  the  noncatalytic  wall  condition  Eq.  (IV-1) 
is  the  applicable  heat  transfer  expression  for  cases  both  with 
and  without  mass  injection.  Furthermore,  a  noncatalytic  wall 
condition  must  be  imposed  on  a  chemically  reacting  flow  in  order 
to  permit  meaningful  comparison  with  nonreacting  flow  heat- 
transfer  results.  This  restriction  has  obvious  importance  in 
interpreting  the  numerical  results  presented  in  this  report 
and  comparison  with  experimental  wind-tunnel  data. in  a. chemi¬ 
cally  reacting  flow. 
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is  abstract  The  effects  on  tjje  pressure  and  heat- transfer  distributions  of 
injecting  helium  and  argon  into  an  air  laminar  boundary  layer  over  a 
spherically  blunted  7.5-deg  half-angle  cone  at  =  13.2  and  Re^/in.  = 
1520  were  considered.  Calculations  are  presented  for  a  binary  perfect 
gas  boundary  layer  including  transverse  curvature  and  displacement  ef¬ 
fects,  an  ionizing  air  nonequilibrium  boundary  layer  without  higher- 
order  effects,  and  a  nonequilibrium  stagnation  point  thin  viscous  shock 
layer  model.  The  effects  of  inviscid  chemical  equilibrium  and  non¬ 
equilibrium  external  conditions  and  wall  catalyticity  were  investigated. 
Neglecting  all  higher-order  effects,  comparison  of  binary  gas  and  ioniz¬ 
ing  air  boundary- layer  results  showed  reasonable  agreement.  Strong  ef¬ 
fects  of  inviscid  chemical  nonequilibrium  on  stagnation  point  heat  trans¬ 
fer  were  found.  Arguments  are  given  for  using  inviscid  nonequilibrium 
stagnation  conditions  rather  than  the  usual  chemical  equilibrium  con¬ 
ditions  for  nonequilibrium  boundary- layer  calculations. 
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